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Введение. Предлагается метод решения задачи о бесконеч- 
ной пластине, лежащей на упругом основании. На пластину 
действует периодическая нагрузка в виде силы, перемеща- 
ющейся по произвольной замкнутой траектории. Цель ис- 
следования — разработка численного метода решения задач 
теории упругости для тел, находящихся под действием по- 
движной нагрузки. 

Материалы и методы. Учитывая периодичность рассмат- 
риваемой нагрузки, она раскладывается в ряд Фурье на вре- 
менном отрезке, длина которого равна периоду нагрузки. 
Решение исходной задачи строится посредством суперпози- 
ции решений задач, соответствующих нагрузке, задаваемой 
слагаемыми описанного выше ряда Фурье. Окончательное 
решение задачи представляется в виде отрезка ряда. Каждое 
слагаемое при этом соответствует решению задачи о воз- 
действии на бесконечную пластину нагрузки, распределен- 
ной по замкнутой кривой (траектории движения силы). Для 
нахождения этих решений используется фундаментальное 
решение уравнения колебания бесконечной пластины, ле- 
жащей на упругом основании. 

Результаты исследования. Предложен новый метод реше- 
ния задач теории упругости для тел с нагрузкой, движущей- 
ся по замкнутой траектории произвольной формы. Решена 
задача о бесконечной плоскости, по которой с постоянной 
скоростью движется сосредоточенная сила. Определено, что 
траектория движения представляет собой гладкую замкну- 
тую кривую, состоящую из дуг окружностей. Рассмотрен 
характер изменения перемещений и напряжений вблизи 
движущейся силы. Изучено распространение энергии упру- 
гих волн. С этой целью выполнено вычисление координат 
вектора Умова — Пойтинга. Исследовано влияние скорости 
движения силы на длину вектора Умова — Пойтинга. 
Обсуждение и заключения. Метод применим и при рассмот- 
рении более сложных объектов (плиты сложной формы, 
слоистые плиты, вязкоупругие плиты). Его преимущество 
— экономичность, так как для построения решения исполь- 
зуются уже известные решения задач. Окончательное реше- 
ние выражается в удобном виде — как сумма криволиней- 
ных интегралов. Полученные результаты могут быть ис- 
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пользованы в процессе проектирования дорог. Изучение 
распространения энергии упругих волн от движущихся 
транспортных средств позволит оценить воздействие ука- 
занных волн на строения, расположенные вблизи дороги. С 
учетом данных о характере изменения перемещений и 
напряжений оценивается износ дорожного покрытия. 


Ключевые слова: бесконечная пластина, движущаяся 
нагрузка, произвольная замкнутая траектория, энергия 
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Введение. Исследование динамических явлений, обусловленных действием подвижной нагрузки, пред- 
ставляет собой актуальную задачу, имеющую прикладное значение (например, при решении вопросов развития 
транспорта). В работах, посвященных указанной проблеме, рассматривались различные задачи с подвижной 
нагрузкой. В частности, показано, как на полуплоскость или полупространство (упругое изотропное, трансвер- 
сальное изотропное, вязкоупругое) действует нагрузка, движущаяся по бесконечной прямой линии с постоян- 
ной скоростью. В этом случае при решении задачи вводится подвижная система координат, связанная с дви- 
жущейся силой, что позволяет исключить время из числа независимых переменных [1-5]. Часть работ посвя- 
щена изучению действия на бесконечную пластину или полосу (упругую или вязкоупругую) равномерно дви- 
жущейся по прямолинейной траектории нагрузки. При этом используется тот же метод исключения временной 
переменной или рассматривается квазистатическая постановка задачи [6-12]. Наибольший интерес представ- 
ляют задачи, в которых длина траектории движения нагрузки конечна, а сама траектория представляет собой 
кривую линию. Часто этом случае используется конечноэлементное моделирование движущейся нагрузки [1 |- 
13]. В ряде работ при решении подобных задач используются вариационные методы (в частности метод Рэлея 
— Ритца) [14—16] или разновидность метода Галеркина, позволяющая сводить задачу к обыкновенным диффе- 
ренциальным уравнениям. При этом рассматриваются разнообразные объекты приложения подвижной нагруз- 
ки (плиты, слоистые плиты, вязкоупругие плиты, полупространства — как изотропные, так и анизотроп- 
ные) [17—19]. В данной статье представлен метод, который развивает идеи, изложенные в работах [20—23]. 

Материалы и методы. Рассмотрим бесконечную пластину, лежащую на упругом винклеровском осно- 
вании, которая находится под действием нормально приложенной силы, перемещающейся по замкнутой траек- 
тории. 

Задача сводится к интегрированию уравнения движения пластины, лежащей на упругом винклеровском 
основании [14]: 


ты 
АИ +0 Е =—, (1) 

р 

ЕЙ? 
где И’ — прогиб пластины; ДО = 12—) ; Е — модуль Юнга; у — коэффициент Пуассона; й — толщина пла- 
—У 
> рй К 

стины; с “= г ; р — плотность материала; К = г ; К — коэффициент податливости упругого основания; Р — 


сосредоточенная сила, движущаяся по замкнутой кривой 7 с постоянной скоростью а. 
Введем координату 5, отсчитываемую от некоторой фиксированной точки кривой у. Тогда сила Р, дви- 
жущаяся вдоль кривой '/ со скоростью а, будет описываться соотношением Р=Р(5-—аё). Функция Р(5-—ай) 


рр 
периодическая по 1, с периодом Т = — ‚ где Ё — длина кривой у. 
а 


Решение. Рассмотрим установившийся режим. Разложим функцию Р(5-—а!) в ряд Фурье по перемен- 


ь [ 
ной {на отрезке | —-—; — |. В этом случае коэффициенты разложения предстанут виде: 


2а 2а 
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Г 
т 
акт, 
с: = | Ро -ае 1. 
_Г 
а 
Сделав замену переменной интегрирования в интеграле 5—аё = 2 , получим: 
Г, 
я 215 225 
7 е У. е У 
с. = | Р(еде а, 
[@ 


-7 


А, в 
где 4, — коэффициенты разложения в ряд Фурье на отрезке |= | : функции Р(2). 
Тогда подвижную нагрузку можно представить в виде ряда Фурье: 


1 — 2ип(5—а!) 
Ре-а) = У Че й | 
а 
К=—осо 


Учитывая линейность задачи, ее решение можно представить в виде: 


Я 


со 
И = А.И’, 5 (2) 
К=Ь—со 
где И’, — прогибы пластины, обусловленные действием вертикальной нагрузки, распределение которой вдоль 
К 


25 —Р ати 
кривой у описывается функцией е ^/, изменяющейся во времени по закону е Я 


Для определения И’, воспользуемся фундаментальным решением уравнения (1), которое соответствует 


кап 
т 





Р=5(х-ж)5(0?- уе" ‚ где о} = 


Используя принцип предельного поглощения и традиционные методы построения решений дифферен- 
д 
циальных уравнений, можно получить фундаментальное решение уравнения (1), которое при А > имеет 
б 
ВИД: 


и Ух, У) = [№ (%)- Ко (а>К) › 


1 2 | 
где К = |= } +(У- у у № = ^ 6, ; © = ве = о. ‚ Ко(2) — функция Макдональда. 
йе 


2 


© 
При А < —— решение имеет вид: 
с 


1 т 
= — НН (ув)- К, (78 


2 
— -К, Н ры (УК) — функция Ганкеля. 





- 21п(5—аг) 
Тогда И’, = фи (х, у, ху (5), У (5))е и. 45. Используя известные формулы теории тонких пластин и 


у 
полученную по приведенным выше соотношениям формулу, определяющую прогиб И’ (2), можно вычислить 


перемещения и, и, и напряжения с., о, и ох, в любой точке пластины. 
При больших А приходится вычислять интеграл от быстро осциллирующих функций. Для этого исполь- 
зовалась квадратурная формула, основанная на замене кубическим сплайном слабо осциллирующей части 


2Кт5 
подынтегральной функции, а сильно осциллирующий множитель е *^” рассматривался как весовая функ- 


ция [15]. 
Результаты исследования. Итак, бесконечная пластина лежит на упругом винклеровском основании. 
На нее действует нормальная сила, перемещающаяся по траектории, изображенной на рис 1. 
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Рис. 1. Траектория движения нагрузки 


Расчеты проводились при следующих исходных данных: й = 0,25 м; с=221 м/с; Е= 232469-10 Н/м’; 
у= 0,36. На рис. 2 приведены результаты расчетов, соответствующие К = 1,864 м“ иа= 125м/с. 
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Рис. 2. Изменение перемещений и напряжений 


Данные графики описывают изменения перемещений и напряжений в системе координат, связанной с 
точкой приложения движущейся силы Р. Ось # направлена по касательной к траектории движения {‚ аосьи — 


по внешней нормали к области, ограниченной траекторией п (см. рис. 1) при 2 = р 
При этом вектор перемещений и тензор напряжений соответственно представлялись в виде 
И =Ш-1 п.п +’. , 5=51.И +51. пт+5т-(т+пГ), 


где А — нормаль к пластине. 
На рис. 2 показано изменение вдоль оси { компонент вектора перемещений ИТ, ИИ, Ци! ‚ тензора напря- 
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жений 51, ИЕ, УЕ и изменение этих величин вдоль оси и — Ии,т,Иипи и эт, пи, “ти. 

Согласно расчетам, характер поведения указанных величин во всех точках траектории остается неиз- 
менным. Расчеты показали также, что при изменении скорости движения нагрузки а в пределах от 0 до 125 м/с 
компоненты вектора перемещений и тензора напряжений возрастали незначительно (на 3—4%). 211 
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Для изучения распространения энергии упругих волн вычислялись компоненты вектора Умова — Пой- 
тинга е; =—5, 0; (5, — компоненты тензора напряжений, и; — координаты вектора перемещений), направ- 


ление которого указывает направление распространения энергии, а длина описывает количество энергии, пере- 


носимое через единицу поверхности, перпендикулярной к направлению данного вектора в единицу времени. 
На рис. 3. представлено распространение энергии упругих волн около движущейся сосредоточенной 


силы (положение силы обозначено красной звездочкой). 





Рис. 3. Распространение энергии упругих волн 


Обсуждение и заключения. Анализ полученных результатов показывает, что при указанных выше 
пределах изменения скорости длина вектора Умова — Пойтинга практически пропорциональна скорости движе- 
ния нагрузки и качественная картина распространения энергии вблизи силы слабо изменяется в процессе дви- 
жения. Характер поведения рассчитанных выше перемещений и напряжений во всех точках траектории остает- 
ся неизменным, а их значения слабо зависят от скорости движения нагрузки а при изменении этой скорости в 
пределах от 0 до 125 м/с. 

Применение предложенного метода к бесконечной пластине, лежащей на упругом основании, не ис- 
черпывает его возможности. Он может быть использован при рассмотрении более сложных объектов (плиты 
сложной формы, слоистые плиты, вязкоупругие плиты). Рассмотренный метод отличается от упомянутых выше 
большей экономичностью, так как использует для построения решения уже известные решения задач. Оконча- 
тельное решение выражается в удобном виде — как сумма криволинейных интегралов. 

Полученные результаты могут быть использованы в процессе проектирования дорог. Изучение распро- 
странения энергии упругих волн от движущихся транспортных средств позволит оценить воздействие указан- 
ных волн на строения, расположенные вблизи дороги. С учетом данных о характере изменения перемещений и 
напряжений оценивается износ дорожного покрытия. 
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